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Abstract. The subfactor approach to modular invariants gives insight into 
the fusion rule structure of the modular invariants. 



1. Introduction and Background 

We are going to use the inter-relations between quantum and classical dynam- 
ical systems and to use tools from non-commutative geometry or non-commutative 
operator algebras to understand phenomena in classical statistical mechanics. 

A good starting point is the two dimensional Ising model and its C*-treatment 
[H I23L I25| . The classical model is set in the configuration space V = {±}^ 
of distributions of -I- and — on the two dimensional lattice Z^, and the nearest 
neighbour Hamiltonian 

a.J3 nn 

where the sum is over nearest neighbours a and (3 in Z^, and a — (cr") in V. 
The transfer matrix method allows us to study the classical model set in C{P) = 
(^22 and its equilibrium states, characterized by say the Dobrushin-Lanford- 
Ruelle equations or a variational principle by a quantum system of noncommuting 
observables A — ^^M2 in one dimension with dynamics at = )T^'* and 

associated equilibrium states or more precisely its ground states. Associated to an 
equilibrium state /i at inverse temperature (3 for the Ising model is a ground state 
(Pf^ on the Pauli algebra A and to each local observable F in C{V), a quantum 
observable Fp depending only on the temperature such that we can describe the 
classical correlation values in terms of quantum ones: /i(i^) = (^^(i^/j). The key 
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Figure 1. Yang-Baxter equation 



element in this reduction is to identify the classical partition Z function in the 
denominator of the classical expectation value 

(1.1) ii[F) = ^('^) exp(-/3i?(a))/ exp(-/3H(a)) 

in quantum terms. We write for horizontal periodic boundary conditions: 
Z = E.exp(-/3i7(a)) 



= E 



= Y T^iiiJ^i^i^ ■ ■ ■ ^?iv«i = trace T...T = trace = trace e"^^. 

a 

The partition function of a square of size M by N decomposes as the exponential 
of a sum factorising into a product of the partition functions of columns summed 
over column configurations and hence can be understood as a matrix product, 
computed in the trace as we have imposed periodic boundary conditions. Taking 
the scaling limit at criticality, we obtain a field theory, and for a torus model with 
periodic boundary conditions in addition vertically, we have 
(1.2) Z(t) = tr(q(^«-^/24)q(£o-c/24)) 

Here Lq and Lq are commuting Hamiltonians, usually part of commuting Virasoro 
algebras as for the Ising model here with associated central charge c, or multiplier 
of the projective representation of the vector fields Im = —z"^^^d/dz on the 
circle. Much of this structure can be understood or is present in the statistical 
mechanical model itself, and one of our aims is to lift structures or understanding of 
the conformal field theories at criticality back to the original statistical mechanical 
models. For example, the central charge c itself can be seen from the asymptotics of 
the partition function Z: when 1 << iV << M, then Z « exp{-NMf + MiTc/N6), 
where / is the free energy log{Z)/MN after N,M ^ oo US]. To generahse the 
Ising model we need some integrability of solutions of the Yang-Baxter equation as 
in Fig. n for local Boltzmann weights. Such solutions are naturally provided by 
Hecke algebras or quantum groups associated to SU(n)k particularly at roots of 
unity - so that SU{2)2, the Ising model is the first in a double series of examples. 
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The Boltzmann weights he in 

((g)M„)^'^("^' =7r(Hecke). 

The justification of the term SU{n) models is as follows. By Weyl duality, the 
representation of the permutation group in A/„ is the fixed point algebra of 
the product action of SU{n). Deforming this, there is a representation of the 
Hecke algebra in Af„ whose commutant is a representation of a deformation of 
SU{n), the quantum group SU{n)q |31) . The braid relations are then precisely the 
YBE at criticality. The representations which appear here are labelled by Young 
tableaux of at most n — 1 rows but a further constraint implying rationality or 
finiteness of the representation labels occur at roots of unity q — e'rV("+fc) when 
only labels with at most k columns appear. The SU{2)k models are constructed 
by distributing edges of the Dynkin diagram A^+i on square lattice. The first non 
trivial example here is described by A2, the two extreme vertices representing our 
previous symbols {±} and the internal symbol is a dummy variable. The labels thus 
have two meanings coming from the representation theory of SU{n) of symmetric 
group (and their deformations) and which point of view we want to emphasise 
may depend on whether we are more interested in the statistical mechanics or the 
conformal field theory picture: 



statistical mechanics 


conformal field theory 


Hecke algebras 


SU{N) loop groups 


A symmetric or braided 
group representation 


A positive energy 
representation of SU{N) 


III bimodules 
Jones- Wenzl subfactors 


nil sectors 
Jones- Wassermann subfactors 


Ocneanu paragroup 
and connections 


Longo Q-systems and 
canonical endomorphism 



2. Modular Invariants 

The link between the two frameworks are deformed Weyl duality, Popa's classi- 
fication of type III inclusion in terms of those of type II and of course the partition 
function Z . Decomposing Eq. H1.2|l according to our underlying loop group sym- 
metries we have: 

— ^ A,/^ 

where Xx{t) = tr/f^ exp(27ririo) is the character in the positive energy represen- 
tation A and q = e^'^" . Modular invariance under reparameterisation of the torus 

then becomes Z{t) = Z {{ar + b) / (ct + d)) for (^yj in SL{2]1) which is generated 
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by 5 = Q^) and T = (J J). Since 5*^(2; Z) acts linearly on characters, there is a 
representation of SL{2\ "L) hy S ^ S = [^a,^], and T — [Ia,^], where 

Xa(-1/t)=V 5'a,/.XmW> Xa(t+1) = V Ta.mXmM- 
Then modular invariance is 

(2.1) ZS ^ SZ.,ZT = TZ. 
Since the coefficients Z\^ appeared as multiplicities, 

(2.2) Za^G 0,1,2,..., 

and usually through uniqueness of the vacuum we have 

(2.3) Zoo = 1. 

A modular invariant is then a matrix Z = [Za^] satisfying the above constraints 
Eq. H2.ll) . Eq. (|2.2|) and Eq. H2.3|l . although we will sec that much interesting struc- 
ture is uncovered if we relax the normalization condition Zoo = 1- Usually as in 
the SU{n)k case, the modular data is such that S is symmetric and ^oa > 0. 
In this case P |28, 

(2.4) < Za,^ < dxd^ 

where d\ ~ Sq,x/Soo, and thus there are at most finitely many normalized modular 
invariants. 

There always exists the trivial modular invariant Z = 1. For SU(2) there is 
the celebrated A-D-E classification [14], e.g. SU (2)i6 have three modular invariants 
labelled by the Dynkin diagrams A17, Dio, E7 with corresponding invariants Z^n = 
1, Zdio is an orbifold obtained by folding the Aiy-diagram which has a fixed point, 
and Ze7 is a twist of the orbifold invariant. These are the normalized modular 
invariants, and if one relaxes the normalization condition. Zoo = !> to consider for 
example products of normalized modular invariants then an interesting fusion rule 
algebra appears as in Table 





ZAiv 


ZUlo 


ZEr 


ZAi7 


ZAiv 


ZUlo 


Ze-, 


Zdio 


Zdio 


2Zdio 


2Ze, 






2Ze, 


Ze>io + Zet 



Table 1. Fusion of SU {2)iq modular invariants 



It is then natural to ask whether the (unnormalized) modular invariants form a 
fusion rule algebra generated by the normalized modular invariants and get a better 
understanding of these fusion rule algebras. The subfactor approach can assist us 
in this [^1^ . 

3. Moore Seiberg dilation 

Here we will have a type III factor N , endowed with a system A of endomor- 
phisms which are taken to be braided and so yield a modular data. As we have 
seen we always have the trivial modular invariant 

(3.1) ^ = EaI^^I'' 
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where we now interprete x\ ^-s a formal character trq^"^'^/^^ even when the Hamil- 
tonian Lq may not exist. More generally we may introduce twists or permutations 
of the fusion rules preserving S and T and the vacuum so that we should also 
consider 

(3-2) Z = ^^x\X*i){xy 

In some sense every modular invariant can be dilated or brought to this form. 
Suppose we can extend the system in the following sense. We have a subfactor 
N C M with a system B of endomorphisms of M. One can emphasise the extension 
aspect of moving from endomorphisms of N to those of M (a- induction j34L E]). 
One needs the system A to be braided and the sector of Tl, the dual canonical 
endomorphism to lie in T,{A), the set of finite integral sums of endomorphism of A. 

Under these conditions we emphasise the restriction aspect of this extension. 
The trivial or twisted modular invariants Eq. H3.2I) for the ,B-system can restricted 
to A, written formally in terms of characters as Xt — "^xbr.xXXi with branching 
coefficients &t,a- Then restricting the diagonal modular invariant 

(3.3) = \Xr\', 

to the original system we have: 



(3.4) Y.\xr\' = J2\Y,K,xx\ 

reB tEB XeA 

with mass matrix 



(3.5) Zx.f_, = br.xbr,fj,. 

These invariants are called type I (or more precisely the inclusion N G M describing 
this modular invariant is type I) and are necessarily symmetric Zx^ — Z^^x- In the 
presence of non-trivial twist d of the B system, we have type II invariants 



(3.6) Zx.,_i — br,xb^(T),t^- 

which have symmetric vacuum coupling Zqx ~ Zxo- 

A type III phenomena accurs, due to some underlying heterotic structure which 
results in needing different labellings B+ and B- on left and right extended systems 
A C B±. In the subfactor framework, this is found in two intermediat e subfactors 
N C M± C M, where M± carry systems of endomorphisms B±. Not only can this 
situation be found with modular invariants with non-symmetric vacuum coupling 
as in the orthogonal at loop groups low levels jH], but also with quantum doubles 
of finite groups '26 . Turning to the extension point of view N C M, when A C 
End(-/V) is non-degenerately braided and the dual canonical endomorphism 9 lies in 
Ti(A), we form to be the systems of endomorphisms on M from the irreducible 
components of the sectors of {a"^ : AG A}. These generate the full system C = 
C+ V C~ , not necessarily even commutative, but B = = C+ n C~ is not only 
commutative but non-degenerately braided |S]. 
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4. Izumi quantum Eq model 

Non-degenerately braided systems and the corresponding modular data can be 
obtained in the operator algebra setting in at least two ways. One is to take the 
Jones- Wassermann loop group examples in algebraic quantum field theory analysed 
by Wasserman ?43 and students e.g. Laredo- Toledano 141) . The other is to take 
the quantum double of a system of endomorphisms - which may not even be com- 
mutative let alone braided or may be braided but the braiding has some degeneracy. 
Here we focus on one such example the quantum double of the even Eq system. 
(As explained in |13l Section 5] , this is more natural than taking the double of the 
entire Eq system as the Longo-Rehren construction then gives the double of the 
Eq subfactor). The Eq system can be realised through a conformal embedding of 
SU{2)io C 50(5)1 HHISIII]. The modular data obtained from the Longo-Rehren 
inclusion of the even system of Eq is explicitly written down as follows in '3U', Page 
648] where d = 1 + VS, X = 2 + (f ,i = The S and T matrices are 



/ 1 


1 


1 + d 


1 + d 


2 + d 


d 


d 


d 


d 


d X 


1 


1 


1 + d 


1 + d 


-2 - d 


d 


d 


-d 


-d 


-d 


1 + d 


1 + d 


1 


1 


2 + d 


-d 


-d 


-d 


-d 


-d 


1 + d 


1 + d 


1 


1 


-2 - d 


-d 


-d 


d 


d 


d 


2 + d 


-2 - d 


2 + d 


-2 - d 




















d 


d 


-d 


-d 





d 


d 


-d 


~d 


2d 


d 


d 


-d 


-d 





d 


d 


-d 


-d 


-2d 


d 


-d 


-d 


d 





-d 


d 


— 2i — di 


2i + di 





d 


-d 


-d 


d 





-d 


d 


2i + di 


— 2i — di 


l) 


\ d 


-d 


-d 


d 





2d 


-2d 










g(l--l 


1. -1. 1. 


Tri/3 47r 

e , e 


/3 ^5„i/6 


e ' , 


-i). 











4.1. Verlinde matrices of the quantum Eq modeL We computed numer- 
ically the Verlinde matrices of the quantum double of iJg : Nq is the identity matrix, 
and the others can be written as the following quadratic expressions 

^1 = XoXl +XlXS +X2X3 +X3X2 + IX4p + X5X6 +X6X5 +X7X8 +X8X7 + IX9p, 
^^2 = X0X2 +X2X0 +X1X3 +X3X1 + 1X2 + X3p + (X2 +X3)(X5 +X6)* 

+ (X5 + X6)(X2 + X3)* + 2|X4P + X4(X7 + X8 + Xo)* + (X7 + X8 + X9)X4 
+X5X6 + X6X5 + X7(X8 + Xg)* + (X8 + X9)X7 + X8X9 + X9X8: 
^^3 = X0X3 +X3X0 +X1X2 +X2X1 + 1X2 +X3p + (X2 +X3)(X5 +X6)* 

+ (X5 + X6)(X2 + X3)* + 2|X4P + X4(X7 + X8 + X9)* + (X7 + X8 + X9)X4 

+ |X5|' + Ixel' + |X7|' + |X8|' + X9(X7 + X8)* + (X7 + X8)x;, 
Ni = (X0 + X1+2X2 + 2X3)X1 + X4(X0 + X1 + 2X2 + 2X3)*+X4(X5+X6)* 

+ (X5 + X6)X4 + (X2 + X3 + X5 + X6)(X7 + X8 + X9)* 
+ (X7 + X8 + X9)(X2 + X3 + X5 + Xe)*, 
^^5 = (XO + X3 +X5)X5 +X5(X0 + X3)* + 1X2 +X3p 

+ (X4 + X7 + X8 + X9)X4 + X4(X7 + X8 + Xg)* + (Xl + X2 + X6)X6 

+X6(xi + X2 + Xe)* + X7X8 + X8X7 + Ixgl^, 
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^6 = XoXe +X6X0 + X1X5 +X5Xl + 1X2 +X3p +X2X5 +X5X2 +X3X6 

+X6X3 + (X4 + XT + X8 + X9)X4 + X4(X7 + XS + Xg)* + X5X6 

+X6X5 + IX7|' + |X8|' + |X9|', 
^7 = X0X7 + X8X0 + XlX8+X7Xl+X2(X4 + X8+X9)* + (X4 + X7 + X9)X2 

+X3(X4 + X7 + Xq)* + (X4 + X8 + X9)X3 + (X5 + X6)X4 + X4(X5 + Xe)* 

+X5X8 + X7X5 + X6X7 + X8X6: 
^8 = XOXS +X7X0 + X1X7 + X8X1 +X2(X4 + X7 + X9)* 

+ (X4 + X8 + X9)X2 + (X4 + X7 + X9)X3 + X3(X4 + X8 + X9)* 

+ (X5 + X6)X4 + X4(X5 + Xe)* + X5X7 + X7X6 + XeXs + X8X5> 

^9 = X0X9 +X9X0 + (X2 +X3)(X4 + X7 + X8)* + (X4 + X7 + X8)(X2 +X3)* 

+ (X5 + X6)(X4 + X9)* + (X4 + X9)(X5 + Xe)* + X1X9 + X9Xl- 

All the iVi's are symmetric apart from N7 and Ng that are transpose of each 
other. The first two matrices Nq and A^i are permutation matrices as the primary 
fields and 1 are the simple currents of our present modular data. The quantum 
Eq modular invariant Z2 is the charge conjugation invariant, i.e. [Z2]^ ^ = S\,p,- 
The Frobenius-Schur indicator FSx = 1 for all A G ^ except for A = 7, 8 where it 
is zero. 

4.2. Modular invariants of the quantum subfactor(s). The dimen- 
sion of the commutant {S, T}' is four and spanned by the modular invariants which 
are exactly four and were computed numerically using 10, Eq. (1.3)]. We obtain 
Zi = Nq, and the others are given by the following quadratic expressions: 

Z2 = |Xo|' + |Xl|' + IX2|' + |X3|' + |X4|' + |X5|' + |Xe|' + |X9|'+X7X8 + X8X;, 
Z3 = |X0+X2|' + |X1+X3|' + 2|X4|', 

Zi = |xo+X2+X4p- 

4.3. Subfactors for quantum Eq modular invariants. Here we consider 
the Longo-Rehren N C M inclusion of the even system of the Eg Dynkin diagram, 
and let A be the system of endomorphisms that yield the above modular data. 
By [261 Lemma 3.8] (see also |3Dl) we can choose an endomorphism Ai from the 
sector [Ai] such that A^ = 1 since T\^^Xi = 1- Therefore [9] = [Aq] © [Ai] is a dual 
sector of the braided subfactor TV C x Z2 . The irreducible decomposition of the 
N-M sectors is as follows: [tAo] = [tAi], [iA2] = [tAs], [LX4] = [ai] © [02], [/.As] = 
[tAg], [tAy] = [tAs], and [tAg] = [03] © [04]. Hence gives rise to a trace 8 modular 
invariant. So ZNcNyiZ2 — Z2 the permutation invariant. Second the LIi{Ee) dual 
canonical sector is [6] = [Xq] © [A2] © [A4]. The irreducible decomposition of the 
N-M sectors are as follows: [iAo], [iAi], [iA2] = [iAo] © [lX^], [lX^] = [lXi] © [lX^], 
[6A4] = [lXq] © [tAi] © [iAs], [lAs] = [lXq] = [iXj] = [lXs] = [tAg]. Therefore LR(£;6) 
produces a trace 3 modular invariant. Hence Zi^Yi{Ee) = ^4- By ^13, Sect. 5] the 
notation of Fig. 1 on page 28, [(0,0)] © [(6,0)] is the [Aq] © [A2] in our current 
notation. In |13[ Example 5.1] instead of taking the three even vertices of Eg 
one can take the two extreme endpoints (which are even and discard the internal 
point) following Izumi's Galois correspondence |29l Proposition 2.4]. In this way, 
we get an intermediate subfactor N <Z P <Z M with inclusion map Lp such that 
[9] = [Ao] © [A2] is the sector of the dual canonical endomorphism of A C P. Then 
by the Verlinde matrices we find that: [tpAo], [tpAi], [tpA2] = [ipAo] © [ipAs], 
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[ipXi] = [ipAg] [a] ® [6], [tpAs] — [ipXe], [lpXt] — [ipAs], [tpXg], which gives trace 
6 modular invariant. Thus the subfactor N C P produces Z3. Moreover, chiral 
locality holds true here. We thus have proven the following. 

Theorem 4.1. All the modular invariants of the quantum double of the Eq 
system are realised by braided subfactors. 

The fusion rules between the quantum Eq double (sufferable) modular invari- 
ants are as in Table El It is worth mentioning that as the modular invariants form a 
basis of the commutant {S, T}' the coefficients exist and are uniquely determined, 
but due to a fusion rule algebra phenomenon those coefficients are non-negative 
integers [261 Corollary 3.6] (see also |2H). 





Zi 


Z2 


Z3 


Zi 


Zi 


Zi 


Z2 


Z3 


Zi 


Z2 


Z2 


Zi 


Z3 


Zi 




Z3 


Z3 


2Z3 


2Zi 


Zi 


Zi 


Zi 


2Z3 


3Zi 



Table 2. Fusion ZaZl of modular invariants 



[ao] [ai] 




[ay] [ag] [ag] 



Figure 2. Zi, fusion graph of [a^] 

4.4. Full systems of the quantum Eg modular invariants. The global 
index of our system is oj — 8{l + d + d"^) ^ 89.5692. 

Case Z\. For the trivial modular invariant Z\ we have ^ ~ C*' ~ C. We display 
the fusion graph of \ct^ in Fig. |21 

Case Z2. For the permutation invariant Z2, the full system C is obtained by 
permuting the sectors of those in the case Z\. 

Case Z3. For the sufferable quantum Eq modular invariant Z3 we have by |12j : 

Lu±^ 8(1 + d + d^)/{2 + d)^ 18.9282, luq = 4. 
Computing using the Verlinde matrices we find: (ao,Q!o) — {af,a^) — {af,af) 

= (a^>"^> = ("7 ."t") = ("s'.^s) = 1. (a^>a^> = ("f = 2, {a^,af) = 3, 
(ao, Oif) = 6o^i + 52,t, {af,af) = 6i^i + (53,i {af,af) = (a^, a^) = {at,a^) = 

{at,af) = (a^,a6> = 1' ("^,"7) = ("t^as) = (aJ^ag) = 1, {af,af) = 
(a^,af) — (af^ag) — (ag ,ag) = 1, and the others vanish. We then conclude 
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Figure 3. Z3, fusion graphs of [af] and [af] where Z| = 2Z3 



that [ao] and [af] are irreducible, and [af] — [ao] [af], [af] — [af] [af], 
[af] = [af] [af^'^] [af^\ [af] = [af], [af] = [af] = [af]. Using also 
the matrix Z=Z3 and noting that [a^^*''] — [aj^*''] for i = 1,2 we conclude that 
C" ~ {ao, ai, a'"^\a''^''}, and — U {af,af}. Moreover since 

{afa^,af) = nI^Zi^^ = 3 = (aja^', 04 aj") 

we have the sectors of C° equal Z2 x Z2 with [q!4 ■*] [04*''] = [id]. We need to find 
four more irreducible sectors in the full system C because #C = Tr(ZZ*) = 12. We 
compute and find that 

{afa^,afa^) = ^ 7V« gTV^'^Z^,, = 2 

and the dimensions of the intertwiner spaces between afa^ and any sector in 
both chiral systems vanish. Hence [a^a^] decomposes into two new irreducible 
sectors [(a^a^)^^)] and [(a^Q;^)^^)]. The other two sectors are similarly obtained 
from the decomposition of afa^ further noting that {af a'^ , af a^) = 0. We 
denote these new irreducible sectors by [{afa^ and [{afa^Y'^^. We have the 
fusion graphs of both [af] and [a^] displayed on the LHS of Fig.|3| We use straight 
lines for the fusion graph of [af] and dashed lines for that of [a^]. On the RHS of 
Fig.|31the fusion graph of [af] and [a^] is displayed. We display the fusion graph of 
[af] and [a^] in Fig.^ When 9 is the local dual canonical endomorphism Aq A2, 
let 7 denote the canonical sector of C M. Since {afa'^,j) — (a^,a^) = 
1 we conclude from 7. Corollary 3.19] that [7] = [ao] [(afa^)^''-'']. The full 
system C decomposes as two sheets. The first is the chiral system C+ and the 
second sheet comprises the irreducible components of C^a^ in accordance with the 
decomposition of Z| — 2Z3 and [7] = [ao] [(^5^0:5^)^^^] and j26. Corollary 3.6] 
and in]. 
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Figure 4. Z3, fusion graph of [af ] where Z| = 2Z3 



Case Z4. The global indices of the chiral systems and C° are by il2) : 

w± = 8(l + d + d2)/(4 + 2d) = 2 + (i2 - 9.4641, = t^±/t^ = 1 

respectively. There are 9 = Ti{Z4Zl) M-M irreducible sectors in the full com- 
mutative chiral system C. Note that C*^ = {ao} since ujq = 1. We next com- 



The fusion rules are as the well known fusion rules of the even vertices of the 



graph Ef, m = [«o] ® [af] ® 2[af], [af][af] = [af], [af][af] = [ao]. 



We find that C = C+ x C and remark that d^± = 1 + V3. We display the fusion 
the all M-M system together with the fusion graphs of both [a^] and [af] in the 
LHS Fig.[S| In this figure, we use straight fines for the fusion graph of [af ] whereas 
dashed lines for that of [a^]. We also encircled the C^-chiral sectors with small 
circles and with larger circles those for the C~-chiral system. The RHS of Fig. [3 
we display the fusion graphs of [af ]. If 9 is the local dual canonical endomorphism 
Ao ® A2 ® A4, let 7 be the corresponding canonical endomorphism of C M. As an 
application of [3 Corollary 3.19], we can compute that [7] = [ao]® [af aj"]©[af a^]. 
In accordance with Z| = 3^4 and j26[ Corollary 3.6] the full system C decomposes 
as tfiree sheets C"*", C"*'aj~ and C^a^ . 
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